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ABSTRACT 

This document discusses the application of physical-optics diffraction 
theory to a deployable dual-reflector geometry selected for study by the 
Jet Propulsion Laboratory, California Institute of Technology, Pasadena, 
California. The methods employed herein are not restricted to the 
"Conical-Gregorian" antenna, but apply in a general way to dual and even 
multiple reflector systems. Complex -vector wave methods are used in 
the Fresnel and Fraunhofer regions of the reflectors. Field amplitude, 
phase, polarization data, and time-average Poynting vectors are obtained 
via an IBM 360/91 digital computer. Focal-region characteristics are 
plotted with the aid of a CALCOMP plotter. Comparison between the GSFC 
Huygens-wavelet approach, JPL measurements, and JPL computer results 
based on the near-field spherical wave expansion method are made where- 
ever possible. 
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Symbol 


Meaning 


F 

a 

max 

°0 

a 

c 

z i 

f 

R 

C 

d, D 

3 

r, 6, 

P, 0, <t> 
N 

E, H 
x', y\ z’ 
x, y, z 
n. 

i 

V 

co 

¥ 

P, e, <? 

ds 

k 


the i th surface when tracing rays or waves 
focal length of the paraboloid 
maximum radius 
minimum radius^ 

inclination angle of a parabolic arc 

inverse slope of a cone 

displacement of cone apex from origin 

wavelength 

frequency 

Rayleigh far-field criterion 

diameter 

feed function 

observer coordinates radius, theta, phi 

coordinates used with feed functions 

exponent in feed function (pos., zero, or neg.) 

electric and magnetic vector fields 

observer position in Cartesian coordinates 

surface coordinates 

vector unit normal to surface 

del operator 

angular frequency 

solution to wave equation 

constitutive parameters: magnetic permeability, inductive 
capacity, electric conductivity 

differential area 

wave number 
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time-average Poynting vector 
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AN ANALYSIS OF A DUAL-REFLECTOR ANTENNA SYSTEM 
USING PHYSICAL OPTICS AND DIGITAL COMPUTERS 

INTRODUCTION 

The purpose of this document is (1) to provide a written record of a GSFC 
analysis of dual and multiple reflector systems, (2) discuss some aspects of 
the Conical-Gregorian arrangement, and (3) make a comparison between the 
GSFC Huygens -wavelet approach, JPL measurements, and JPL near-field 
spherical-wave expansion methods. Data accumulated by JPL for the Conical- 
Gregorian antenna is not reproduced here. 1 

The geometry of the system is shown as Fig. 1, together with the input 
parameters. It should be noted that a part of the conical main reflector is 
not illuminated under ray optics. An illumination function 9 = cos N 0 , N = 12.3750 
was used in the GSFC studies to achieve an edge taper of -18.0 db on the sub- 
reflector edge. Under ray optics the -18.0 db taper maps to the interior of the 
cone rather than the edge. A parametric representation is used for both surfaces, 
and can be found in the Appendix with other details of the computation. 

FORMULATION OF THE SCATTERED FIELDS 

Consider the formulation for the scattered fields (under transmission) and 
allow that every point of the conical main reflector is deep inside of the Fresnel 
region of the subreflector. Assume, further, that the observer of the fields of 
the main reflector may lie in either the Fresnel or Fraunhofer regions of the 
latter. 

The fields due to the main reflector are computed via a formulation which 
is equivalent to Kirchhoff-Kottler and Franz formulations. 2 


^ef 1, page 146 
2 Ref. 2, Chapter 8, pp. 460-470 
Ref. 3, Chapter 5, pp. 158-162 
Ref. 4 pp. 141-144 
Ref. 5, pp. 114-115 
Ref. 6, pp. 500-506 
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E (x 
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H (x\ y', z') = 
where 

and 


— f ( n x HL) x Vfds, 
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¥ = 


» _ i k r 


V 


ip = -( j k+ l)wi r 


taking l r in a local context on the reflecting surface. 

Since only the magnetic field H 2 is required on y 2 , the integral 


H 2 (x', y', z') = - J_ f (HjxH^xVfds 

77 J^i 

over y 1 is sufficient to obtain the illumination of the conical reflector. 

Some studies on the subsystem were performed in the Fresnel and Fraunhofer 
regions. A separate program was used in this instance to evaluate 


E ( 
and 


x\ y\ z') = _ — _ f [ n xH.)^ VNs-j w/x — C n ! * H, 

] V e 4 TT I 1 1 4 77 I 

J v. 


) V^dS 


H (x', y', z') 


-M 

4 Tt I 

J -V 


(Hj x Hj) x V W d s 


It is noted that every value H^ on y 2 is obtained by summing the Huygens 
wavelets over all of y (an m:l mapping) which is distinct from ray optics (a 1:1 
mapping). Scattered fields under reception are obtained in the same manner as 
those under transmission, however, the conical reflector is now y x and the 
subreflector is now y 2 . With this interchange, the preceding integrals for 
E(x\ y\ z') and H(x’, y\ z') in the focal region can be formed by the same logical 
process. 
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THE SUBSYSTEM 


On occasion it is helpful to study a subsystem to effect design changes, obtain 
better insight with regard to the electromagnetic fields and wavefronts, or to 
verify satisfactory performance and avoid programming errors in a multi- 
reflector system. The results of the GSFC far-field computations for the 
subreflector are given in Fig. 2 and 3. These agree quite closely with the 
measured and calculated data {not reproduced here). The central spike at 6 =0° 
is observed and is about 17 db below the peak at 9 =50° (the inclination angle a 
of the parabolic arc). The single cut provided by the JPL article shows a central 
spike about 18 db below the peak at approximately 6 - 50°, reading the graphical 
results. Feed displacement (p e ) equals zero here. 

Although good agreement was obtained when comparing the far-field sub- 
reflector patterns of two independent simulations and a measurement, these 
patterns do not impinge on the conical main reflector. For this reason, JPL 
utilizes the near-field spherical wave expansion method. In the GSFC approach 
the Huygens wavelet approach is retained, but consideration is given to the near- 
field geometry in summing the wavelets. Actually, it would appear that the time- 
average Pounting vector provides a more comprehensive measure of the electro- 
magnetic effects in a Fresnel region than either the electric or magnetic fields. 

The basis of this viewpoint is the fact that power flow is affected when fields 
are asynchronous in time and non-orthogonal in space. The time -average Poynting 
vector incorporates all of these effects, to whatever degree they may be present, 
since 


<P> = i-Re ExH* 

2 

is sensitive to spatial orientation of the fields and the relative phase between 
electric and magnetic fields. This product is also sensitive to the magnitude of 
each of the fields, which may be especially important if the electric and magnetic 
fields in a near-field or intermediate near-field region are not related simply, 
as in a far-field region. 

Plots of the time-average Poynting vector passing through the aperture plane 
of the conical main reflector are given as Fig. 4 and Fig. 5. It can be seen that 
the power density of the central spike is of the same order of magnitude as the 
intensity at 6 =50°, which is distinct from the 18 db difference shown by the 
far-field pattern discussed previously. Solid lines in Fig. 4 and Fig. 5 indicate 
the intersection of the geometrical bounds with the aperture plane. 
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Detailed studies of the subsystem were made by exploring areas of the 
incident ray bundle from the subreflector in the vicinity of the conical main 
reflector. These were Fresnel region studies of field amplitude, phase, and 
the time-average Poynting vectors for both scanned and unscanned transmission 
modes. The plots are highly detailed and not presented here. In general they 
established the following. The near -field analysis agrees quite closely with the 
ray-optics predictions in most respects - wavefronts are orthogonal to rays 
from the subreflector, and time-average Poynting vector plots show power-flow 
along the rays. It is noted, however, that the diffraction analysis shows details 
not provided by ray-optics. Rays carry with them no accurate measure of 
intensity, whereas the Poynting vectors were significantly reduced near the 
edges of ray bundles, etc. The power flow for both scanned and unscanned cases, 
as well as the wave fronts, were exceptionally orderly or homogeneous. 

In contrast to this, the electric fields due to the subreflector showed 
interesting characteristics. Constant-intensity plots yielded a lamellar-type set 
of contours in the central portion of the ray bundles. This can be attributed to 
the prime-feed weighting of the subreflector illumination distribution, modified 
by the effects of phase stationarity in the Fresnel region. The computed divergence 
of electric field intensity via diffraction was found to agree closely with the law 
or rule of divergence that can be inferred from energy conservation principles. 

Since the rays due to the subreflector pass orthogonally through the slant 
surface of a truncated imaginary cone, and since total far-field radiated energy 
is invariant, it follows that 


E a [it (2 r tan a - s cos a ] 1/2 

V P P 

When range r from the subreflector is large, the above reduces to 

E a JL 
/7 

which is the same as the law of divergence for cylindrical waves. In the preceding, 
(s) is the width of the ray bundle in Fig. 1, and the slant height of the imaginary 
cone. As before, a p is the inclination angle of the parabolic arc used to generate 
the subreflector surface. The expression has significance for values of r such 
that the power flows through a cone or' a truncated cone. 
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THE DUAL REFLECTOR SYSTEM 


The backscattered radiation patterns of the dual-reflector system 
were computed using the previously cited parameters to determine the agreement 
with JPL measured data and the near-field spherical-wave expansion method. 

The GSFC approach depends entirely on Huygens wavelets. Previous discussion 
pointed out that the illumination on the subreflector, obtained here by a directive 
point-source feed function is a one-to-one (1:1) mapping. The illuminatjon on 
the main reflector is by contrast, a many-to-one (m:l) mapping over y . Fresnel 
or Fraunhofer fields due to y 2 are likewise obtained as a many-to-one mapping 
(n:l) taken over r 2 . This should be contrasted with ray-optics analyses. 

Three sets of integration sampling intervals were chosen for the main and 
subreflector surfaces. These were LI 2 = 3.0, LI t = 1.0; LI 2 = 1.5, LI X = 1.0; 

LI 2 = 1.5, LI t = 0.5 respectively, in wavelengths. An analogue to Shannon's 
sampling theorem indicates that curvilinear quadrilaterals on antenna surfaces 
should not exceed 0.5 \ on an edge, however, it has been found advantageous to 
explore intervals both larger and smaller to obtain a satisfactory compromise 
or trade-off between accuracy and economy. Planar "cuts" at = 0°, 45°, and 
90° were taken for each set of LI specified, above, to establish the stability of 
the solutions with respect to LI. It was found that cpu time was approximately 
3/4 minute, 1 minute, and 3 minutes for the sampling criteria employed here. 
These cpu times should be assessed after recalling that significant deviations 
among the patterns occurred in the vicinity of the third side lobe and beyond. 

Figure 6 corresponds to the cut 4> = 0.0°, obtained with LI 2 = 1.5, LIj = 0.5, 
and is an E -plane cut showing both the magnitude of the electric field | E 0 | , and 
the associated phase 'P E0 . It was found that cross-polarization components for 
principal-plane cuts were 80 db to 100 db below beam peaks and are, therefore, 
not presented here. In the diagonal or <P = 45° cuts cross-polarization was about 
25 db below beam peaks, and the form of the pattern was the characteristic vee- 
shaped cross-polarization pattern that is obtained with single-reflector systems 
(paraboloids, for example). H-plane and diagonal cuts were very similar to the 
E-plane cut, and are not presented here. 

An appreciation of the effect of sampling can be obtained from Fig. 7, which 
also shows JPL measured and calculated results superimposed on GSFC data 
for the 0 = 0° pattern cut. Side-lobe levels are presented as a function of 
integration sampling intervals (LI) up to, and including, the fifth sidelobe. It is 
noted that the convergence to the solution is oscillatory and not monotonic with 
respect to sampling (LI). JPL measured data is double-valued for a given 
sidelobe in some instances and this is probably attributable to surface fabrication 
errors, system asymmetries, and slight pattern-range variations. In any event, 
the largest measured deviation is 3 db, occurring at a level of 26 db below the 
beam-maximum level. 
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A very brief study pertaining to the stationarity of the phase of the waves 
arriving at the main reflector aperture plane was made using LI 2 = 3.0 and 
LI j = 1.0. It should be recalled that the JPL conical main reflector, Fig. 1, ex- 
tended beyond the geometrical interior bound from cr Q2 = 1.40 ft to o- = 1.17 
to enhance the system directive gain. By ray optics there is no access to the 
shaded portion of the conical main reflector (i.e. no illumination from the sub- 
reflector), but all wavelets from y 1 have access to every point on y 2 • The 
following Table sums-up the findings when LI 2 = 3.0 and LIj = 1.0. 

Inner Radius of Conical Electric Field \E 9 \ 

Main Reflector (ft.) (db scale) 


2.00 -76.84 

1.75 -75.38 

1.50 -74.63 

1.40 o- 02 = cr m j (JPL value) -74.77 

1.17 -74.56 

1.00 -74.93 

0.75 -74.84 

0.50 -75.05 

0.00 -74.94 


The large increase in area on y 2 , exceeding the JPL design area by about 
23 percent when o- Q2 = 0.0, results in a decrease in on-axis intensity of about 
0.38 db due to a lack of stationarity. Of the values of cr Q2 used, none lead to an 
intensity exceeding that obtained with the JPL value % 2 = 1.17 ft. (approx.). 

The sharp reduction in gain for o q2 > 1.40 ft is due to the removal of cophased 
or stationary aperture distribution. 


A brief study of the beam scanning characteristics of the Conical- 
Gregorian antenna indicated that beam squint-angle (0 S ) varied almost linearly 
from 0.0° to 0.84° as feed displacement (x € ) in the xy feed plane varied from 
0.0 \ to 1.5 A. Comma-lobe levels about 7.2 db below beam peaks were observed, 
for a series of computations carried out at the three sets of sampling criteria 
employed earlier for unscanned beams, at the limit of scan 0 = 0.84°. 

S 

The next step, after beam-scanning of individual beams, was to use four 
feeds in a conventional amplitude -sensing monopulse configuration. Feeds were 
displaced A/4 in both x and y directions for z = 0, the feed plane, in an initial 
computation to recover the monopulse sum and difference patterns of this dual- 
reflector system. The resulting sum and difference patterns of Fig. 8 and 9, 
taken at 4> = 0.0°, are representative of the results. Other cuts, obtained at 
4> = 45.0° and =90.0° for both telemetry and error channels are not reproduced 
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here. Sum channel sidelobes as high as 9.0 decibels below the axial peak value 
were noted. Cross-polarization in the 4> = 45.0° cuts reached 26.0 db below 
axial peak level for the sum channel. Error channel cross -polarization levels 
were about 40.0 db below the axial peak level of the sum channel. 

Since the dual-reflector system studied here has a magnification factor 
similar to that of a conventional Cassegrain, the feeds were displaced ^/2 in both 
x and y directions for z = 0. Results were very similar to the case of ^-/4 feed 
displacement, although sum-channel gain was reduced by about 2.3 db. The 
error channel slope was not noticeably affected. It appears that a better selection 
of monopulse feeds could be made by means of a careful examination of the feed 
directivity for a given geometry, increasing the number of feeds to produce a 
larger "effective" area (increasing the number of feeds to 16 for example), and 
by a mapping of the Airy disc and ring structure under reception to determine 
the size of the focal region. The mapping of the focal region in two transverse 
planes was carried out for the Conical- Gregorian antenna and is discussed in 
the last section of this report. 

Numerous other studies were either made in detail, or planned, as they were 
within the capability of the program. Some of these are "squinted" monopulse, 
variation of feed directivity, and studies of system sensitivity with respect to the 
axial and lateral placement of the feeds and the reflector elements. These are 
not presented here as they are straightforward applications of the program and 
are not of general interest. 


FOCAL-REGION MAPPING 

Focal-region mapping 1 was restricted to the case of axial plane-wave recep- 
tion. The simulation approximated the latter by means of an isotropic source 
(in amplitude) polarized in the "vertical" sense, P(I) = 1.0 polarization moment, 
and situated 100 miles from the origin of coordinates (the feed point of the 
system). Under reception, the plane wave of the simulation impinges on the 
conical main reflector (7 t ), and is converged to the subreflector (y 2 ) where 
the final convergence is to the focal region. See Fig. 10. Ordinarily the Airy 
disc and ring structure is of interest, as is the depth of field of the system. 

These can be found by searching for, or computing, fields in the xy and yz planes, 
respectively, for the polarization selected here. It is anticipated that the Airy 
disc, whose radius is given by 2 


’Ref. 7, pp. 777-783 
Ref. 8, pp. 935-943 
2 Ref. 9, p. 283 
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Figure 10, Ray-trace for reception 
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for parabolic reflectors, should be broadened due to the magnification factor of 
the Conical-Gregorian system. Likewise the depth of field should be extended. 

The "Airy disc" structure was obtained on the assumption that convergence 
would be in the feed-plane of the system. Only the electric field E x was plotted 
to illustrate the disc and ring structure of the focal region using polarization 
moment P(I) = 1.0 for the incident plane-wave. The results are shown as Fig. 11 
and Fig. 12 with nulls and a few quantitative values retained. These contours 
were obtained from a field of values at a resolution of k / 10 with LIj = 3.0, LI 2 = 
1.0 and LIj = 1.5, LI 2 = 0.5. Inadequate sampling is depicted by Fig. 11. A third 
computer run, not illustrated, at a sampling interval LI X = 1.0 for the conical 
main reflector and LI 2 = 0.5 for the JPL subreflector, established the fact that 
the results of Fig. 12 were stable with respect to the integration interval. 

In order to present a graphic illustration of the behavior of the focal-region 
for the Conical Gregorian antenna in the transverse or yz plane cut, the 
electric field E x , its associated phase 'i' Ex , and the time-average Poynting vector 
<P> are displayed as Fig. 13 to 15. The region mapped is of the same extent 
and has the same scale, inwavelengths, for the preceding quantities. It is noted 
that E z = 0, and E x is all of the main polarization; the cross polarization E is 
vanishingly small. Also, the Poynting vector is composed of P y and P z com- 
ponents only, P x being omitted (intentionally) in the yz cut. In the original plot, 
which has been photographically reduced here for publication, the maximum com- 
puted Poynting vector was limited to L q = 3/8-inch on the plot and the minimum 
value resulting from the calculations was represented as a point or dot. That is, 


P - P . 

L = min 

P - P • 

max min 


where all power values are taken in decibels. A brief analysis of isophote, wave- 
front, and power-density plots follows. 


Figure 10 also illustrates the region selected for study, above, in addition 
to tracing the incoming rays through the dual-reflector system. It can be seen 
that all of the rays converge at a point focus (F) according to the theory of ray 
optics. In the event that there is no detector at (F), there is a divergent cone of 
rays in addition to the convergent cone, as shown. It is noted that these rays 
cross the z-axis at F, the apex of the cones. It is interesting to compare Fig. 10 
with Fig. 13 through Fig. 15. and note similarities and differences between ray 
optics and diffraction theory for the Conical Gregorian system. 
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Figure 13 shows the constant amplitude contours (isophotes) of the electric 
field E x , and provides considerable insight relative to the focussing characteris- 
tics of the system. It can be seen that the null-to-null depth of field at this par- 
ticular frequency, is about 8 wavelengths. The null-to-null width of field is about 
3 wavelengths. (The latter can be obtained from Fig. 12, Fig. 13, and Fig. 15.) 
Even though ray-optics indicates symmetry about the xy plane, diffraction com- 
putations show considerable asymmetry at 8.4 GHz. The focal region is dis- 
tended due to the magnification factor of the system, but is otherwise similar 
to that of paraboloids. 1 

Figure 14 illustrates the constant phase contours (wavefronts) for the same 
region depicted by Fig. 13. It can be seen that the perturbed regions of this 
wavefront plot correspond to the regions where low field intensities (nulls) were 
observed in Fig. 13. It is possible for the wavefronts of different value (phase) 
to join at a null field point without contradiction since phase is undefined when 
amplitude is zero in complex-variable analysis. The wavefront plot should also 
be compared with Fig. 15, the Poynting vector plot, since the time-average power 
flow is everywhere orthogonal to the wavefronts. It appears that an increase in 
phase velocity of about 7 percent is present in the focal region. 

Figure 15 shows that the time-average Poynting vectors, <P> = 1/2 R e E x H*, 
obtained via diffraction theory are very different from ray-optics. The geometric 
bounds shown in this figure are those of the cones of Fig. 10. Power flow tends 
to conform generally to the prediction of ray optics, but with notable exceptions. 

For example, the Poynting vectors are bent, or diffracted to pass through the 
Airy disc and rings without crossing the'system axis. Null regions appear in 
the field, sharp umbral and penumbra regions are everywhere absent, and a 
vortex-like counterflow of power is observed in the zero-energy isophote regions. 

The scaling of the Poynting vectors is such that any noticeable reduction in 
the length of a vector represents a significant change in the power density at 
that point. For a dynamic range of 60 db, a reduction of the maximum length of 
Poynting vector L 0 to 5/6 L 0 , 4/6 L 0 , 3/6 L y corresponds to a reduction in 
power density of 10 db, 20 db, and 30 db respectively. A careful examination 
of Fig. 15 shows that very low power densities exist in the umbral region with 
the exception of those vectors passing through the central portion of the Airy 
disc. Due to the reduction of Figs. 13 through Fig. 15 from joined Calcomp 
sheets measuring about 6 feet by 3 feet to the 8-inch by 10-inch document page, 
it is impossible to extract accurate quantitative data from these reproduced 
figures. The data field values from which the contours were plotted have been 
deleted everywhere except around the borders and a few other strategic locations. 
Photographic reduction would have rendered them unreadable. 

’Ref. 10, pp. 439-449 
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An enlarged or "magnified" illustration for the region of low intensity near 
the point (x, y, z) = (3 A., 0, - X./2) is given as Fig. 16. Here the Poynting vectors, 
isophotes, and wavefronts are presented in a single sheet to assist the reader in 
superimposing these features of the electromagnetic fields. Figure 13 through 
Fig. 15 resolved the field into A/10 increments. Figure 16 resolves the field into 
A /so increments, and shows an unusual flow of power around a weak isophote that 
appears to be at least 55 db below the largest field value observed on the system 
axis at z = - V2. The phase singularity is seen to coincide with the "null" in 
intensity. 

Further "magnification" of the "null" region was attempted by resolving the 
field into A/200 increments. Figure 17 once again shows Poynting vectors, 
isophotes, and wavefronts. A definite vortex-like counterflow can now be seen 
for the Poynting vectors. An even weaker isophote, about 85 db below the largest 
field value can be identified, and the wavefront representation continues to suggest 
a singularity. The persistence of the counterflow and phase singularity through 
VlO, A/50, A /200 resolution tends to clarify the fundamental behavior of the 
electromagnetic characteristics in the vicinity of a null. The results also indi- 
cate that the limits of accuracy have not been exceeded with the computer since 
a high degree of order can be observed in the plots at A /200 resolution. 


SUMMARY 

The results obtained in this investigation by means of the GSFC diffraction 
program agreed closely with measured and computed JPL results for the "Conical 
Gregorian" configuration. Since JPL Fresnel-region computations were predicated 
on the use of the spherical-wave expansion method to obtain the illumination 
distribution on the conical main reflector, it is reasonable to infer the equivalence 
of the GSFC Huygens -wavelet approach to the expansion method. Although no 
data are presently available for comparing the computed focal-region fields 
(isophotes), wavefronts, and time-average Poynting vectors, the results appear 
valid in terms of the ray-optics predictions. The formulation used herein is, 
furthermore, the same as that employed with single reflectors (paraboloids) to 
obtain focal-region information. 1 

Subsequent to the investigation of the Conical Gregorian dual-reflector 
system, at 8.4 GHz, Ku band 60-foot Cassegrain systems (approaching 1000 A 
by 1000 A) and three -reflector systems have been studied by employing slight 
variations of the method described herein. It has been found advisable to study 

’Ref. 11 
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subsystems in detail in the Fresnel region, establish sampling stability (LI), and 
treat the scattered subsystem fields as prime-feed data. Interpolation is often 
admissible between a modest number of pattern "cuts" over the subsystem. 

This approach should be contrasted to the method described in this report, where 
each charge and current value on a surface is computed by integrating over some 
other entire surface without any "a priori" knowledge of stationarity on the latter. 
Symmetry is not assumed with the methods employed in this report. 

Present plans call for extensive application of the existing program to ground- 
based and spacecraft antennas. A number of new subroutines, accommodating 
unique project geometries and even composite surfaces, are being annexed to the 
main program. Long-term plans include simulation for electrically-small scatterers 
(d < 1A), artificial dielectrics, and antenna synthesis. This document is an 
interim report. 
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APPENDIX A 

CONICAL MAIN REFLECTOR GEOMETRY 


x = o' sin ^ 
y = - cr cos £ 
z = c a + Zj 


Tangents to Surface 


_ dp — d p — 7 7 

p r-— v > 


y + k z 


Unit Normals to Surface 

— _ r x _ '°cr x '°5 _ (- c sin £ i + c cos £ j + 1 k) 

|^x^| (EG - F 2 ) 1/2 ~~ JO + c 2 ) 1 / 2 !' : 


Differential Area 

d S = (E G - F 2 ) 1/2 d ct d £ = | ~p a x p J dcrd£ 

E = x 2 + y 2 + z 2 = 1 + c: 2 

F = x^ x^ + y^ y^ + z^ z^ = 0 (orthogonal net) 
G = x 2 + y^ 2 + z 2 = ct 2 

d S = cr (c 2 + 1) 1/2 d a d £ 
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APPENDIX B 

SUBREFLECTOR GEOMETRY 

A parabolic arc 

z" = x" 2 /4 F 

under a rotation transformation 


T x : z" 

= z' 

cos 

a - 

x' 

sin 

a 

x" 

= z' 

s in 

a + 

x' 

cos 

a 


about its vertex becomes 

(z' cos a - x' sin a)/4 F - (z' sin a + x' cos a) = 0- 
A subsequent translation transformation 

T 2 : z' = z + F cos a 
x' = x - F sin a 

which causes the focal point to coincide with the origin of coordinates yields 

z 2 sin 2 a + z (2 x sinacos a. - 4 F cos a) + (x 2 cos 2 a + 4 F X sin a - 4 F 2 ) = 0- 
The quadratic in z then leads to 

= Z= C^F cos a- 2x sin acos a) ± [(2x sinacos a-4Fcos a) 2 -4 sin 2 a(x 2 cos 2 a+4Fxs ina-4F 2 )] 1 2 

z sin 2 a 

The singularity in the last equation implies the trivial case of rotation for a 
parabolic arc. 
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The parabolic arc becomes imaginary when (x) exceeds 


(2 x sin a cos a - 4 F cos a) 2 - 4 sin 2 a (x 2 cos 2 a + 4 F x sin a - 4 F 2 ) = 0 

or when (x) exceeds 


x = F/sin a (a ^ 0) 


The corresponding value of (z) is 

z = Fcos a/sin 2 a (a ^ 0) 


Differentiating the solution (z) for the translated and rotated arc with respect 
to (o'), having let (+ x) - (o'), gives 


z - cot® [2(2o-s inacosa-4Fcosa)(2sinacosa)-(4sin 2 a)(2o r cos 2 a+4Fsina)] 

4s in 2 a [( 2<vs in a cosa-4Fcosa) 2 -4s in 2 a(cr 2 cos 2 a+4o-Fs i na-4F 2 )] 1/2 

which reduces to 


z — — co t a ( 

( 7 


sin a (1 sin a 

F 


1/2 


The two singularities in the equation, above, have been identified in the preceding 
discussion. 

@ a = 0, tt, . . . , it n 

implies the trivial rotation, and 


S 2 @ CT = 


s i n a 


implies the imaginary domain for the parabolic arc. 
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Surface 


x = cr sin £ 
y = - a eos £ 

Z = (a) + Zj 


Tangents to Surface 

— dp— d p _ ■ ~ - - 

p <. = rr = tp p ° 1 * + 1 y + k * 

Unit Normals to Surface 

n = ^ ^ ^ x ^ - z a sin £ 1 + z a cos £ j + 1 k 

|p CT xpJ (EG - F 2 ) 1/2 (i + z^l/2 

Differential Area 

d S = (EG-F 2 ) 172 dad £ = |p«, x p, | dad £ 

E = x z + y 2 + z 2 = i + z 2 

o * cr cr cr 


F = + y a + z a z^ = 0 (orthogonal net) 

G = + y^ + z 2 = a 2 


dS = (E G) 1/2 dad 5 = 

[»♦ 

- co t a e ^ 



L a . \l/2 

sin a 1 sinal 


L 

L l F J 


2 - 1/2 


adcrd(. 
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